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Abstract. The optimal strategy for a microscopic swimmer to migrate across a linear 
shear flow is discussed. The two cases, in which the swimmer is located at large 
distance, and in the proximity of a solid wall, are taken into account. It is shown 
that migration can be achieved by means of a combination of sailing through the 
flow and swimming, where the swimming strokes are induced by the external flow 
without need of internal energy sources or external drives. The structural dynamics 
required for the swimmer to move in the desired direction is discussed and two 
simple models, based respectively on the presence of an elastic structure, and on an 
orientation dependent friction, to control the deformations induced by the external 
flow, are analyzed. In all cases, the deformation sequence is a generalization of the 
tank-treading motion regimes observed in vesicles in shear flows. Analytic expressions 
for the migration velocity as a function of the deformation pattern and amplitude 
are provided. The effects of thermal fluctuations on propulsion have been discussed 
and the possibility that noise be exploited to overcome the limitations imposed on 
the microswimmer by the scallop theorem have been discussed. 

PACS. 47.15.G Low-Reynolds-number (creeping) flows ~ 87.19.ru Locomotion 



1 Introduction 



There has been recently a resurgence of inter- 
est in low Reynolds number swimming. Partic- 
ular attention has been given to discrete designs 
[1-3], which allow simpler description of the ge- 
ometrical aspects of the problem, and identifica- 
tion of the necessary ingredients for propulsion. 

One of the reasons for renewed interest is 
the progress in mechanical manipulation at the 
microscale, which has allowed the realization of 
the first examples of artificial microscopic swim- 
mers [4-7]. These examples constitute the first 
step towards the construction of "microbots", 
whose application would be wide-spread, e.g. in 
medicine, as microscopic drug carriers in not 
otherwise accessible regions of the human body. 
At the present stage, however, most of such ar- 
tificial swimmers are driven by external fields. 



and the problem of an autonomous power source 
remains under study. 

Several solutions to this problem have been 
proposed. Among them, various methods of rec- 
tification of Brownian motion [8,9], and mechan- 
ical reactions in the swimmer body, induced by 
inhomogeneity in the environment, e.g. presence 
of a chemical gradient [10-12]. 

It has recently been suggested, that a mi- 
croswimmer may extract the energy needed for 
locomotion, out of the velocity gradients in an 
external fiow [13]. Based on a discrete design, 
that is a generalization of the three-bead swim- 
mer of [2, 14, 15] , it was shown that a microswim- 
mer could migrate across a linear shear fiow, by 
a sequence of deformations induced by the ex- 
ternal flow itself. A continuous version of this 
microswimmer has been described in [16], based 
on the analogy of the deformation sequence in 
the discrete case, with the tank-treading motion 
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regime of a vesicle (or of a microcapsule) in a 
linear shear flow [17-19]. 

Through tank-treading, a microscopic object 
such as a vesicle, is able to maintain a fixed 
shape and orientation in an external flow, with 
its surface circulating around its interior, pre- 
cisely as a tank-tread [20]. The existence of a 
preferential shape and orientation for the ob- 
ject, turns out to be one of the main ingredi- 
ents for migration in an external flow. It should 
be mentioned that tank-treading has already a 
long history, as a propulsion candidate for swim- 
ming in quiescent fluids [1, 21]. Other exam- 
ples of migration induced by simultaneous rota- 
tion and deformation of solid objects in viscous 
shear flows, have been analyzed both theoreti- 
cally [23] and experimentally [24] . The design of 
a microswimmer that uses its internal degrees of 
freedom to spin rather than to swim, has been 
proposed in [22]. 

The interesting aspect of external flow aided 
propulsion, is that the migration velocity scales 
linearly in the stroke amplitude. This behavior 
is not surprising: the migration velocity of tank- 
treading vesicles in wall bounded flows, scales 
linearly with the deviation from spherical shape 
[25-28]. In contrast, the velocity of a microswim- 
mer in a quiescent fluid, due to the constrains 
imposed by the scallop theorem [1], is charac- 
terized by quadratic scaling [29]. 

In the present paper, the analysis in [13], 
which focused on the behavior of a discrete mi- 
croswimmer in an infinite domain, will be ex- 
tended to the case of a wall bounded fiow. Par- 
ticular attention will be given to identification 
of the deformation patterns associated with mi- 
gration in different flow conditions, and with en- 
ergy extraction from the flow. A simple semi- 
quantitative analysis of the effect of thermal 
fluctuations will be provided as well, along the 
lines of the approach described in [30] . 

The generation of specific deformation pat- 
terns, requires the presence of a control sys- 
tem, modulating the response of the microswim- 
mer to the external flow (after all, this is what 
characterizes a microswimmer, as compares to, 
say, a simple vesicle). The possibility of con- 
trol through modification of the swimmer struc- 
tural properties will therefore be examined, and 
a simple example of control, by braking on the 
swimmer moving parts, will be described in de- 
tail. 

The paper is organized as follows. In Sec. 2, 
the basic equations of the model are presented 




Fig. 1. Rest configuration of the three-bead swim- 
mer. Small case indicates the reference frame rotat- 
ing solidly with the device. 



and in Sec. 3, the results in [13] are briefly sum- 
marized. In Sec. 4, the analysis is extended to 
the case of a wall bounded flow. In Sec. 5, the 
mechanism of energy extraction from the flow is 
discussed. The structural dynamics of the swim- 
mer is discussed in Sec. 6. In Sec. 7 an hypoth- 
esis of control system to achieve migration is 
discussed. In Sec. 8, the effect of thermal noise 
on passive swimming is analyzed. Section 9 is 
devoted to conclusions. Technical details on the 
swimmer behavior in a wall bounded flow are 
provided in the Appendix. 

2 Swimming strategy in free space 

We analyze the behavior of the simple swimmer 
depicted in Fig. 1. Contrary to [2,15], who con- 
sidered a linear device, the three beads in the 
swimmer under study are located, at rest, at 
the vertices of an equilateral triangle of side R. 
The swimmer is placed in in a linear shear flow 

U(X)=U(0) + (0,aXi,0), (1) 

and wants to migrate along the gradient direc- 
tion Xi . 

We are interested in a situation in which the 
system is able to move by internal deformation 
without the aid of external forces. We assume 
the links between the beads in the trimer to 
be immaterial and the bead radii a to be much 
smaller than R. Including terms up to 0{a/R), 
the equation of motion for the device can be 
written in the form: 

±i=lJ{Xi) + Vi{t)+Fi{t)/a, (2) 

where is the coordinate of the i-th bead, 
Fi{t) is the force on bead i by the rest of the 



OUa: Swimming strategies in low Reynolds number flows 



3 



trimer, a = Qirfia, with fi the fluid viscosity, is 
the Stokes drag, and Ui{t) is the flow perturba- 
tion in Xj generated by the other beads in the 

trimer. 

To lowest order in a/R, the flow perturba- 
tion is obtained replacing the beads by point 
forces in the fluid, with intensity equal to the 
Stokes drag exerted by the beads (stokeslet ap- 
proximation [31]); 

UiW = ET(Xi-X,)F,-; (3) 
T(X) ^"''^ 



\± 
UXI 



4o-L|X 



|X|3 



(4) 



where T(Xi — Xj), for i 7^ j, is the off-diagonal 
part of the so called Oseen tensor [31]. 

Linearity of the shear, Eq. (1), and absence 
of external forces, Yli '^i = 0, imply that the 
trimer center of mass X^" = (Xi -f- X2 -f X3)/3 
would move, if one disregarded the flow per- 
turbation, as a point tracer at X'-'": X°''' = 
U(X<=") + (l/3)2^U,(t). Time averaging the 
deviation with respect to U(X'^'^), we obtain 
the migration velocity 



U™9'^ = (l/3)E(Ui) = (Ui). 



(5) 



The velocity perturbation Ui is the sum of the 
stokeslet fields U(X|Xi,Fi), i = 2,3 generated 
in X = Xi by beads 2 and 3. We assume ex- 
change symmetry among the beads, in the sense 
that trimer configurations at angles (/) + 2n7r/3 
are indistinguishable. This mimicks the station- 
ary tank-treading regime of a vesicle, as the 
trimer arms will extend (or contract) in identical 
way, when reaching a given position in the labo- 
ratory reference frame. Exchange symmetry al- 
lows us to rewrite Eq. (5) in the equivalent form: 

Umisr ^ ([U(X2|Xi, Fi) + U(X3|Xi,Fi)]). 

Focusing on the effect of stokeslet 1 on beads 
2 and 3, rather than on the one of stokelets 2 
and 3 on bead 1, allows to understand geomet- 
rically the swimming strategy of the trimer, as 
illustrated in Fig. 2 (wc consider for the moment 
the case of an unbounded domain). 

As shown in [13] [see Eq. (10) below]_, the 
undeformed trimer will rotate in the field U(X) 
with constant frcqiiency, equal to the fiow vor- 
ticity a/2. Because of this, to lowest order in the 
deformation, U(X|Xi,Fi) will be the stokeslet 
field generated by action of the strain compo- 
nent of U on a (fixed) bead at position Xi . We 
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Fig. 2. The swimming strategy: link 23 is con- 
tracted at = 7r/2 (top-left) and stretched at 
= — 7r/2 (top-right); this causes beads 2 and 3 
to sample the stokeslet field by bead 1 [continuous 
line; see Eqs. (3-4)] differently in the two config- 
urations. The unbalance between the two configu- 
rations, caused by deformation, is responsible for 
migration. Dashed lines indicate the strain compo- 
nents of tJ producing the stokeslet in the two cases. 
Bottom figure: decomposition of U(X) into strain 
and vorticity component. 



see that the the X\ component of the stokeslet 
field generated by bead 1 in response to the 
strain component of U is negative or positive de- 
pending on whcth(;r < < tt or tt < < 27r. 
Due to the approximately constant rotation fre- 
quency, the trimer will spend an equal time at 
the two orientations indicated in figure, and, 
in the case of the undeformed trimer, the con- 
tributions to [/^'sr from the two orientations 
would cancel by symmetry. Deformations, how- 
ever, break this symmetry and beads 2 and 3 
will sample components C/i(X2,3|Xi, Fi) that 
are in general of different amplitude. The end 
result is (/7i(X2,3|Xi, Fi)) 0, which leads to 
migration along X\ . 

Quantitative analysis [see Eq. (14) below] 
shows that the deformation sequence illustrated 
in Fig. 2 will lead to migration of the trimer to 
positive X\ (to the right in figure). 

We notice at this point that migration would 
be achieved also if just one arm in the trimer 
performed the stretching-contraction sequence 
in Fig. 2. [For small deformations, it is possible 
to see that the migration velocity would just 
be reduced by 1/3 with respect to Eq. (5)]. We 
thus see that the presence of an external flow 
makes one degree of freedom sufficient for lo- 
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comotion. We recall that the scallop theorem 
would prevent this, in the case of a microswim- 
mer in a quiescent fluid [1] . Similar "violations" 
of the scallop theorem were observed in [32,33], 
in which case, the role of the external flow was 
played by the perturbation generated in the fluid 
by other swimmers. We shall see in Sec. 8, how 
a similar result could be obtained exploiting the 
presence of thermal fluctuations. 
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3 Migration in free space; 
quantitative theory 

To analyze the deformation dynamics of the trimer 
it is convenient to work in the rotating reference; 
Small case will identify vectors measured in the 
rotating reference frame (see Fig. 1). In the ab- 
sence of rotational diffusion, the motion of the 
trimer will be the sum of a translation and a 
rotation in the plane with rotation fre- 

quency Q = (j), where cf), as indicated in Fig. 1, 
is the angle between the two reference frames. 

Wc arc interested in a regime of small defor- 
mations, and will proceed perturbatively in the 
deformation amplitude. In particular, we will 
have, for the bead position in the rotating frame: 
Xj = x-°'-|-Xj, with x^°' giving the bead positions 
for the undeformed trimer: xf-'/R = {l/VS, 
0, 0), xJ.oVi? = -(l/(2v^), 1/2, 0), x[,°Vi? = (-1 
/(2-\/3), 1/2, 0), and Xj accounting for the de- 
formations. We can express the trimer deforma- 
tions in terms of three arbitrary independent pa- 
rameters Zi (that with two translational and one 
rotational degrees of freedom make up the six- 
dimensional configuration space of the trimer in 
the plane X1X2). With the choice 

Xl = —[y3{z2+Z3),Z^s"' - Z2,0j, 



Fig. 3. Deformations of the trimer corresponding 
to 2i > for 1 = 1,2, 3. In the three cases Zj = for 
j 7^ 



X2 = — - VSZ3, -2zi - Z3, Oj , 
X3 = — - VSz2, 2zi + Z2, oj , 



(6) 



the deformation corresponding to just one Zi be- 
ing non-zero, will describe stretching/ contraction 
of the link opposite to bead i, while bead i 
remains fixed in the plane X1X2 (see Fig. 3). 
Assuming a stationary dynamics for the trimer, 
with exchange symmetry between the beads, we 
can write: 



Zi = "^[An COS ncpi + B„ smn(j)i], 



(7) 



where 0i = (/)2.3 = <t>-^ 27r/3. It will be expe- 
dient in calculating the energy balance for loco- 
motion to consider a hierarchy of deformations 
of increasing order in some small parameter. We 
thus write Zi = z-^' -I- z-^' + . . ., with the nor- 
malized deformation amplitude z — (X]i(4^')^) 
playing the role of expansion parameter for the 
theory. Superscripts will be used to indicate the 
order in z at which a quantity is evaluated. 

Perturbative analysis of Eq. (2) using Eqs. 
(1) and (6) gives the result [13]: 



f (0) _ 

^1 ~ 



aaR ( 1 



a/3 V 
and, to next order 



■ sin 2(/>, 



1 



■ cos 



20). (8) 



fW _ aoRf V3r (1)' 

hi - -^{ — [^2 



f(i) 
'12 



X sin 20] — 

aaR ( 1 
2 12 



1 



^(4^' +^2 



'')cos2(/)|, 
) 

) cos 2^1 , 



(9) 



where primes indicate derivative with respect 
to (j)- Notice that we have disregarded the con- 
tribution from the velocity perturbation Uj to 
the force in Eq. (2). Clearly, fj"' is the reac- 
tion force of the rigid trimer against the external 
flow, and f{^' accounts for deformation. Notice 
that the only contribution strictly qualifying as 
swimming is the one proportional to x^^' in f*^' 
[the ^2,3' terms in Eq. (9); sec also Eq. (7)], the 
remnant being better described as "sailing" . 

The same analysis leading to Eq. (9) gives 
for the rotation frequency [13]: 

f2 = I {2 + [\/3 (z<'' - ) sin 20 + (z<'' 
-h z^'' -2z<'')cos20] -FO(z2)}, (10) 
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and from stationarity of the system, we can re- 
place time averages with angular averages: 



{h) 



— r n{m(i>)d^. (11) 



Expanding J7 in powers of z, we can then write 

the average in the above equation as a sum of 
contributions of increasing order in z: (.) ~ 




(.)<^'+...,with(/)("> = (^a)-i/„'^^ r?< 



~ '+ Fig. 4. The lift on an elongated object in a wall 
((/)) f ((/))d(^.bounded flow. Dashed arrows indicate the strain 



Knowledge of the force in the rotating ref- 
erence frame, allows to write for the migration 
velocity, substituting Eq. (4) into Eq. (5): 



Ijmigr ^ /pj^f^ 



(12) 



where Ti = T(x2 — Xi)-|-T(x3 — Xi) and R is the 
rotation matrix back to the laboratory frame: 



= i?22 — COS0; i?2 



-i?i2 = sin (f){13) 



A simple calculation, including terms up to 0{z 
xa/R) gives for the Oseen tensor Ti: T^^ = 
/3{7/2 - [13zi + 29(z2 + Z3)]/8}, - Tfi ^ 
/3(V3/8)(Z2 - zs) and = /3{5/2 + [2zi - 
31(z2 + Z3)]/}8, with P = 3a/i4aR). Substitu- 
tion into Eq. (12) and exploiting Eqs. (8-10), 
gives the final result, to lowest order in a/R and 



'imgr 



256 



13S<^' 



(14) 



thus t/'"^9''/(ai?) = 0{za/R). Let us consider 
the contribution Sisin^i. We see from Fig. 2, 
that the stokeslet field generated by bead 1 in 
response to the strain component of U, has a 
positive or negative component at beads 2, 3, de- 
pending on whether O<0<7ror7r<0< 27r. 
Migration is produced by the deformation in- 
duced symmetry breaking between the two ori- 
entations. 



4 Migration in the presence of a wall 

A solid wall bounding the flow provides the swim- 
mer with an additional mechanism for migra- 
tion. The flow perturbation by the swimmer will 
be the superposition of what would be produced 
in free space, and a wall correction that can be 
expressed as a sum of images and counterim- 
ages of the freec-space perturbation, generated 
alternatively at the wall and at the surface of 



component of U(X). In the case of a sphere, the 
flow perturbation P and of the image I are reflec- 
tion invariant around the plane perpendicular to the 
flow passing at the sphere center, and the transverse 
drift is zero. In the case of an elongated object, with 
the long axis along the stretching direction of the 
strain, the lobes of the flow perturbation and of the 
image will bo tilted upwards, and the image field 
at the object position will have a net component to 
the left. The opposite will occur for an orientation 
at 7r/2 with respect to the one in figure. 



the beads [31]. In the stokeslet approximation 
of Eqs. (2) and (3-4), only the flrst image has to 
be taken into considerations. 

The migration mechanism is illustrated in 
Fig. 4, in the case of an elongated object: the 
image at the wall of the free-space perturba- 
tion, because of the asymmetry of the configu- 
ration, has a component at the object center, 
that pushes it away from the wall. Now, a rigid 
object, with the exception of the very elongated 
structures described in [34] , will rotate because 
of the flow vorticity. Because of this, such an ob- 
ject will typically alternate between a condition 
of outward and inward drift with respect to the 
wall, and the transverse migration velocity will 
be zero. 

In the case of deformable objects, e.g. vesi- 
cles, a fixed orientation and a non-zero trans- 
verse drift can be achieved by means of tank- 
treading [17]. 

In the case of the triangular trimer, tank 
treading could be mimicked by means of cyclic 
contraction of its arms: the arms contract when 
rotation leads them into the contracting (ex- 
panding) quadrant of the external strain. This 
will produce an overall elongated shape oriented 
along the expanding (contracting) strain direc- 
tion, and lead to a non-zero transverse migra- 
tion velocity. To determine the contribution to 
migration from presence of the wall, it is nec- 
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essary to determine the image of the stokeslet 
field of the beads in the trimer. 

Let us suppose the wall to be located at coor- 
dinate Xi = L with respect to the trimer center 
of mass. Let us denote by U^(X|Xi, Fj) the im- 
age of the stokeslet field U(X|Xi, F;) generated 
in free space by the z-th bead. In Eq. (12), we 
thus have to add a wall contribution: 



As in the case of Eq. (14), we can check that, in 
the absence of deformation, AU^^^^ = 0, and 
that the lowest order contribution with respect 
to z to Eq. (17) is {M^0[x[lfi'^+xi'ifi°^]y°K 

Prom the expressions of Mq,^, and f{°' [see Eqs. 

(8) and (18)], we sec that only harmonics of or- 
der 2 and 4 in z contribute to Z\?7™*^' . Direct 
calculation using Eqs. (6-7) and (8-14), gives in 
fact, to lowest order in z and R/L: 



For small i?/L, we can Taylor expand U^(Xj |Xj, 
Fi) around X^j- = 0. From linearity of low Reyn- 
olds number hydrodynamics, we can write U^(X 
|Xi,Fi) = l(X|Xi)F,;. with I some tensor, and 
the lowest order contribution in R/L will be 
'(0|0)(EjF»> = 0. As regards the first order 
terms, from Xj — 0, we have, identifying 
the three bead contributions with the one from 
bead 1: (E,- X,' • VU^(X|0, Fi))x=o = 0. We 
thus remain with a wall contribution to migra- 
tion: 

^^m^gr ^ g^^^i • V)U^ (0|X, Fl ))x=0 

+ 0{{R/L)'), (15) 

and, from U^(X|Xi,Fi) = l(X|Xi)F„ we ex- 
pect a result in the form AUJ^^^^ = Hap-yiXijiFi^) 
(summation over repeated vector indices is as- 
sumed) . 

In order to determine the coefficients Ha0-y, 
we must know the image field in Eq. (15). The 
image field of a stokeslet induced by a solid wall 
boimding the flow was calculated in [35]. Its 
derivation is outlined for reference in the Ap- 
pendix. We find for the image field derivatives 
entering Eq. (15): 



axi«7/(OlX,Fi)lx=o 
9x,t^/(0|X,Fi)]x=o = 



9aFii 
16^' 



729S<'' 



aaR 
2048L2 • 



(19) 



(16) 



32c7L2 • 

Substituting into Eq. (15), leads to the result: 



rrngr 



27a 
16^ 



(17) 



where Mq,^ = {-RiaRii3 + ^R2aR2i3)- From Eq. 
(13): 



M, 



I cos2(/), 
sin 2^, 



4^4 



COS 2^ 



.(18) 



Notice that the 0{{R/Lf) behavior in Eq. (19) 
is that of the image of a stresslet field at the 
trimer position [31] (the quadrupole field de- 
picted in Fig. 4). As can be seen in Eq. (19), it 
is possible that cancellations arise between the 
n = 2 and n = 4 harmonics, in which case the 
next order in the expansion in R/L for AU"*^^^ 
should be taken into account. Now, the image 
field of a stokeslet at the wall, can be expressed 
as a multipole expansion, whose lowest order 
terms are associated with stokeslets placed at 
the other side of the wall. It is worth mention- 
ing, therefore, that cancellations similar to the 
one discussed above have occur also in the in- 
teraction of two linear swimmers [36]. 

The situation in Eq. (19) corresponds to the 
picture of fixed orientation and migration in- 
duced by tank-treading described at the begin- 
ning of this section. To fix the ideas, consider 
S4 = 0, and focus on the deformation associ- 
ated with zi. The regime i32 < corresponds 
to migration away from the wall. At the same 
time, for B2 < 0, (j) = ±7r/4 will correspond re- 
spectively to contraction and stretching of the 
side 23 of the trimer. In other words, migration 
away from the wall corresponds to the trimer 
maintaining a deformed shape, with long axis 
along the stretching direction of U. This is the 
same behavior of a tank-treading vesicle in a 
wall bounded flow, in the limit of zero viscosity 
of the internal fluid and zero membrane fric- 
tion [25,26]. 



5 Energy balance 

Let us calculate the energy needed to perform 
the swimming actions described in Sees. 2 and 4. 
The average power exerted by the swimmer on 
the fluid is P = 3(xi • fi). This provides a lower 
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bound on the power actually expended by the 
swimmer, which must include the contribution 
by internal friction. 

Let us consider first the case of an ideal trimer. 
The lowest order contribution to the expended 
power, from = 0, is P*^' = (i'^' • f<«))(«). 
Now, in the free-space case of Eq. (14): x^^' = 
x(i)(^/'-ee^ = x/'■^^ where 0/''"" = Bisin<^i + 
B3 sm3(j> and x^^' (z/'^'^'^) is given by the first of 
Eq. (6). Thus, while x^^'^'^ is odd with respect to 
(j), the corresponding force f is even [see Eq. 
(8)] and P*^' = automatically. 

A simple calculation shows that P'^' = 
also for the wall contribution of Eq. (19). Also 
in this case, the mechanism is easy to under- 
stand: focusing on the sequence of contraction 
and stretching of side 23, produced by a defor- 
mation Zi'^'-'' = B2 sin(20) + i34sin(4(/)), wc see 
that the elongation of side 23 for going from 
— 37r/4 to — 7r/4 (energy gained from the fluid) 
is compensated by contraction in passing from 
— 7r/4 to 7r/4 (energy lost to the fluid). In sim- 
ilar way it can be shown that the contribution 
to (x-^''^^ • f<°') and (x™«" • f(«), from Q^^\ in 
the angular average of Eq. (14), is zero. 

The power expended by an ideal swimmer in 
free space, would be therefore 

Pf'""' = 3{xf''''^ ■ ff'""') +0{z^), (20) 

where f-'''^^ = i['\z('''"') [see Eq. (9)], and a 
similar equation is obeyed by the power P""^'' 
that would be expended by the trimer in the 
case of a wall bounded flow (notice that (x-^'"*'® • 

In order for the constrain force f{°' to pro- 
duce work, it is necessary that Xj^' has a compo- 



nent x^"^ = xi(z| 



). with 



= A2 cos 2(j)i. 



From Eqs. (6), (8) and (9), this would corre- 
spond to the trimer extracting from the fluid a 
power 

pextr ^ _^^extr . f(o)^(o) ^ ^^^^ (21) 

To understand the mechanism of power extrac- 
tion, let us focus on the contraction and stretch- 
ing of arm 23 produced by the deformation zf^*''. 
We see that a positive A2, from Eq. (7), corre- 
sponds to link 23 being stretched when it is par- 
allel to the flow direction, and contracted when 
it is perpendicular to it. Energy extraction from 
the flow comes from the fact that the IX2.3I in- 
crease or decrease depending on whether the re- 



spective beads lie in the stretching or contract- 
ing quadrants of the external strain (see Fig. 2). 

In principle, a swimmer could use the mech- 
anism outlined above to extract energy from the 
flow and store it for later use, say, through a sys- 
tem of springs. In alternative, this energy could 
be utilized to compensate the power dissipated 
in swimming, as accounted for in Eq. (20). 

In the case of the ideal trimer described in 
Eq. (20): P'^^*'' = P^''^^ + P""*", which gives 
^extr ^ 0{z^), and A2 = A^f + O(z^) Internal 
friction, however, may contribute to dissipated 
power to 0{z), and a deformation component 
^extr q£ same amplitude as the swimming 
stroke [see Eqs. (14) and (19)] would in this case 
be required. As it will be illustrated in the next 
section, this is going to be a rather natural sit- 
uation, if some kind of elastic structure for the 
trimer is assumed. 



6 Dynamics of the elastic trimer 

We would like to understand the structural dy- 
namics of a trimer undergoing the deformations 
described in the previous sections. 

Let us analyze first the behavior of an elastic 
trimer, in the absence of any internal system 
of control of the device response to the flow. 
Indicating with Vi the angle between the arms 
joining at bead i and with Xij the length of arm 
ij, the potential energy due to bending and to 
stretching can be written in the form 



kbB? 



(22) 



i>j 



where il^i = 17/3+ Aipi and Xij = R+Axij; kbR^ 
is bending elasticity of the joints between the 
trimer arms and Ks is the stretching elasticity of 
the arms. Stretching and bending as a function 
of Zi are obtained from Eq. (6) : 



^2^32 = {-^Zi +Z2 + Z3)R, 

Aiji ^ V3[zi - ^{Z2 + Z3)], (23) 

and cyclic permutations. Substituting into Eq. 
(22), we find the expression for bending and 
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stretching energy: 



AUe 



AUs 



r 2 1 2 I 2 



{ZxZ2 + Z2Z3 + ZzZi)], 



KsR 



■2 r 



33 



{zl + zl + zl) 



+ 12(zi22 + Z2Z3 + 2321) 



(24) 



corresponds to the tank-treading regime with 
long trimer axis along the stretching direction 
of U described in Sec. 4, while A2 > Q corre- 
sponds to the energy transfer from U to the 
trimer dynamics discussed in correspondence of 
Eq. (21). 

The solution to Eq. (27) can be written in 
alternative form as Zi = —\JA^ + B2 sin(2^i -|- 
arctan ^2/52), i-e.: 



Energy balance requires that AUb + AUs + 
AW + AW^n = where AW = Y.^^^ ' 
is the work exerted by the trimer on the fluid 
and AWin is the work against internal friction 
forces. As discussed in the previous section, AW-^^ 
= fj*"' •Z\x,^^' averages to zero in a cycle. From 
Eqs. (6), (8-9) and (11). wc obtain: 



aaR^ 



^^jsin 2(f)i. 



(25) 



Differentiating ZWF*^' + AUb + AUs = with 



respect to Zi 



1, 2, 3, gives the force balance 



equation in the absence of dissipation: 

(6kb + 11ks)2i + {8ks - 3kb){z2 + Z3) 
= sin 201 (26) 

and cyclic permutations for 02,3- Let us assume 
for simplicity that the friction forces acting in 
the trimer are linear in tpi and Xij. Including 
friction leads to an equation in the form 



a 



[(67B + njs)z[ + (87s - 37s) (4 + 4) 



+{6kb + 11ks)^i + i&Ks - 3kb){Z2 + Z3) 
= -^sin2<^i, (27) 

where 7b i?^ and 75 arc bending and stretch- 
ing friction coefficients. Notice that, if 7b, s ~ 
oikb,s, internal friction will produce an 0{z) 
contribution to the dynamics and energy dis- 
sipation in the fluid can be disregarded. 

To solve Eq. (27), we assume a solution in 
the form Zi = A2 cos2(pi + B2 sm2(j)i and obtain, 
after little algebra: (k^2 + ajB2/2) cos2(f)i + 
[kB2 — 07^2/2 -I- 2a(j/9] sin 2(6,; = 0, where k = 
3kb + ks and 7 = 375 + 7s ■ Setting the coef- 
ficients in front of cos20j and sin20i indepen- 
dently equal to zero gives A2 = — 07/ (2^)52 
and B2 = -2ko-q;/(9(k^ -I- q;^7^/4)); in other 
words B2 < and A2 > 0. Notice that i?2 < 



-2a(T 



X sm 



9^k2 + 0^72/4 

• (2(0, ' 



„ arctan^ I ). 
2 2kJJ 



(28) 



With the aid of Figs. 1 and 3 we can under- 
stand the deformation pattern described by Eq. 
(28), and notice the analogy with the behavior 
of a tank-treading vesicle [17] or a microcap- 
sule [19] in a viscous shear flow. In the absence of 
dissipation, the trimer would maintain its long 
axis aligned with the stretching direction of U. 
Adding dissipation would cause the long axis 
to rotate towards the flow, and to get aligned 
with it in the limit a"//K — > 00. No transition 
to a tumbling regime exists. In order for such a 
transition to occur, it would be necessary that 
the trimer rest shape were not that of an equi- 
lateral triangle. Notice that this is the behavior 
of a microcapsule whose rest shape is that of a 
sphere [19]. 

From the analysis in Sec. 4, we see that, in 
the absence of an internal control system, the 
only migration pattern of our trimer, could be 
migration away from a wall bounding the flow. 



7 Swimming through bral(ing 

We have seen that the tank-treading regime of 
Eq. (28), which leads to migration away from 
a wall, is a condition that does not require the 
presence of a particular control system in the 
trimer. Things change if we wish to implement 
the behaviors leading to Eqs. (14) and (19), i.e. 
migration in an unbounded flow and migration 
towards a wall. It has been shown in [13] that 
a simple orientation dependent "braking" sys- 
tem is sufBcient to produce the deformation se- 
quences required for migration. We want to an- 
alyze here the energetics of the system. 

For the sake of simplicity, consider kb — 
Ks = 0, so that the dynamics is dominated by 
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friction, and set 75 = 875/8, so that the sys- 
tem of equations (28) becomes diagonal. This 
is likely not to lead to maximum efficiency in 
terms of speed vs. deformation amplitude, but 
provides an example that is easier perhaps to 
implement experimentally, than variable strength 
springs at the trimer links and joints. 

Under the present hypotheses, Eq. (27) be- 
comes 



4C7 



sin 2(f)i 



(29) 



The two situations leading to drift away from a 
wall, and migration in an unbounded flow would 
require Zi = i?2 sin 2(/)i + . . . with B2 > and 
2, = i?i sin^i + . . ., respectively. 

The first situation could be realized with 
7 = 70 (1 + csm4:(j),)-\ < c < 1. The "brake" 
is released while vertex i has passed the direc- 
tion of maximum stretching, it is acted on en- 
tering the contracting quadrant, and is released 
again after passing the direction of maximum 
contraction. Substituting into Eq. (29) we get 
in fact: 



Zi = 



2a 
970 



[cos 20 j 



1 1 

-I- c{-sm24>i + —sm6(l)i)]. 



(30) 



The second situation could be realized instead 

with 7 = 70(1 + csiiKpi)^^ , \c\ < 1. In this 
case, the brake acts when the vertex is in the 
direction of the flow and is released when it is 
oriented opposite to it (or vice versa, if c has 
opposite sign). Substituting into Eq. (29), we 
get in this case: 



Zi = — ^ [cos 2(/)j + c (sin ( 
970 



isin30i)].(31) 



Notice in both Eqs. (30) and (31), the term 
DC cos20j, that signals energy transfer from the 
fluid to the trimer. 



8 The effect of thermal fluctuations 

Going to sufficiently small scales, thermal fluc- 
tuations will play an increasingly important role 
in the microswimmer dynamics (see e.g. [8, 9, 
30]). 

In the present passive swimming scheme, ther- 
mal noise can act in essentially two ways: 



— Modification of the microswimmer orienta- 
tion with respect to the external flow. 

— Direct contribution to the microswimmer de- 
formation and hence to propulsion. 

The magnitude of the flrst effect depends on the 
ratio of the rotational diffusion time Tf (called in 
[30] flipping time) and the hydrodynamic timescale 
a . (For the sake of simplicity we shall re- 
strict the analysis to the case of a swimmer 
constrained to the plane X1X2). The flipping 
time of the trimer is obtained from the trans- 
lational diffusivity D of an individual bead sus- 
pended in the fluid, from the dimensional re- 
lation Tf ^ /D. The diffusivity of a bead of 
radius a and density p, in a fluid at temperature 
T, can be written in the form 



D ■ 



KT 



(32) 



where K is the Boltzmann constant. Prom here, 
we obtain for the flipping time 



D 



kT 



(33) 



and the condition of negligible rotational diffu- 
sion becomes 

The dimensionless number Pe, that is the Peclct 
number for the bead, parameterizes the relative 
importance of thermal diffusion and hydrody- 
namic transport at scale a. 

Let us pass to analyze the contribution of 
thermal fluctuations to the swimming strokes. 
Let us decompose the deformation in determin- 
istic and fluctuating components: 

6R = 6R + 5r. 

The condition of negligible noise contribution 
to the swimming strokes will be therefore (5r^) 
< {6R^). 

The evolution of the deterministic compo- 
nent 6R is governed by Eqs. (26) and (27). Let 
us consider a situation in which elastic and fric- 
tion force in the trimer contribute at the same 
level to propulsion. This implies a force balance: 



rSR ~ oj^SR R > SR 



(35) 



where F ^ 7/^1, co'^ ^ k/tti and Ta ^ m/a, 
with m ~ pa'^ the bead mass; k and 7 are the 
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(time dependent) elastic constant and friction 
coefficient of the trimer arms; Ta is the Stokes 
time of the beads. Notice that the force balance 
Eq. (35) implies 



Tra ~ - > 1 

z 



(36) 



where Sa = axo <C 1 defines the bead Stokes 
number. 

The evolution of the fluctuating component 
can be described in terms of a stochastic differ- 
ential equation in the form [30]: 

Sr+'^Sr^D'J^ D^^^, (37) 

where ^(t) is a normalized white noise term: 
= S{t - t')- The diffusion coefficient 
Da is determined from the equipartition condi- 
tion 

KiSr"^) ~ m{Sf'^) ~ KT, 

and we can write D/xa ^ DaF ~ t;^^, where Vth 
is the thermal speed of the bead in a fluid at 
temperature T. 

From Eqs. (36) and (37), wc sec that the re- 
laxation time for 5r is the hydrodynamic time 
a^^, and the fluctuation amplitude is {5r^) ^ 
Dar/uj"^. From Eqs. (36) and fluctuating, we 
find the condition of negligible fiuctuating com- 
ponent in the swimming strokes: 



{5E?) 



1 / a 



Pe< 1 



(38) 



A similar approach could be adopted also in the 
absence of an elastic component in the trimer 
dynamics. In this case, the condition of negli- 
gible deformation fluctuations would take the 
form Daa^^ ^ i^R'^), corresponding to the re- 
quirement that diffusion in Sr be small at hy- 
drodynamic timescales. It is possible to see that 
this leads again to the result in Eq. (38) . 

Equations (34) and (38) provide the neces- 
sary conditions for the noise-free analysis in the 
previous section to be valid. We have to ver- 
ify that the resulting migration [Eq. (14)] is not 
overcome by diffusion. 

Over sufficiently small distances, the trimer 
migration will always have a diffusive charac- 
ter. The degree of noisiness of the trimer mi- 
gration could be parameterized in term of the 
ratio of the crossover distance Rr above which 



the trimer trajectory becomes ballistic, and the 
trimer size R. This crossover is easily shown 

to occur at Rc ^ 0/1/™^"^, corresponding to 
a crossover time Tc ^ Rc/lf^^^''' . (Wc are dis- 
regarding the contribution to diffusivity from 
random swimming, which can be shown to be 
appropriate provided the geometric constraints 
z,a/R<^ 1 arc satisfied [37]). From Eq. (14), we 
have t/'"*^'' ~ aaSR/R, and, using Eqs. (32-38), 
we obtain 



f^'^Y p. " ^ Pol/2 



(39) 



and Rc R provided Pe is not too large and 
the conditions in Eqs. (34) and (38) are satis- 
fied. This crossover in space corresponds to a 
crossover in time, which, as described by the re- 
lation 



(40) 



[see Eqs. (34) and (39)], may occur at Tc 

or Tc <C a~^, depending on the relative size of 

the parameters Pe and z. 

Comparing Eqs. (34-39), we see that the ef- 
fect of noise on propulsion can be disregarded, 
if either (or both) Pe and a/SR are very small. 



8.1 Large thermal fluctuations and the 
scallop theorem 

From Eqs. (32) and (34), we sec that the Peclet 
number Pe becomes 0(1) for values of a in the 
range of the micrometer, and shear strengths of 
the order of 1/s. The transition to a diffusivity 
dominated regime is rather sharp due to the in- 
verse cubic dependence of Pe on the bead size 
a. A trimer with submicron beads, would there- 
fore diffuse in the fluid without being able to 
exploit the presence of an external flow; in or- 
der to achieve propulsion, if migration by swim- 
ming remains the strategy, some internal motor 
would be required to overcome bead diffusion. 

It is interesting to sec that thermal noise 
could be exploited to allow the swimmer to pro- 
pel itself with just one degree of freedom avail- 
able. The mechanism turns out to be not partic- 
ularly efficient. Nevertheless, it provides another 
instance of how the constraints of the scallop 
theorem could be bypassed if some additional 
means (in this case rotational diffusion) allows 
the swimmer to change orientation between one 
stroke and its reverse. 
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Fig. 5. Active swimming in a noisy environment. 
Thermal fluctuations cause the trimer to change 
continuously orientation. When bead 1 is pointing 
upward (left) an internal motor stretches link 23; 
the opposite occurs when bead 1 is pointing down- 
ward (right); 52,3 indicate lobes of the stokeslet gen- 
erated by beads 2 and 3. The end result is a velocity 
perturbation Ui at bead 1, with a net component 
pointing downwards independent of orientation. 



The mechanism is described in Fig. 5. An in- 
ternal motor stretches and contracts one link in 
the trimer, when the opposite bead is backward 
or forward oriented with respect to the desired 
direction of motion. 

To fix the ideas, suppose that the fixed frame 
is oriented with the X2 axis pointing upward in 
figure, so that Ui = (0,?7i2,0). The migration 
velocity can be estimated from the velocity per- 
turbation at bead 1 when link 23 is parallel to 
the Xi axis. If ±V is the instantaneous velocity 
of beads 2 or 3 in the deformation, the pertur- 
bation at Xi will be the sum of the two stokeslet 
fields U(Xi 1X2,3, F2,3) = U(Xi 1X2,3, ±(tV). The 
perturbation strength is obtained from Eqs. (2- 
4): 



U12 = 2T2i(Xi3)F3i 



3^3 a 
~R 



V, (41) 



and the sign remains the same inverting the ori- 
entation. (The first index in U12 and labels 
the bead; the two indices in T21 are vectorial). 

As anticipated, this swimming strategy is 
not very efficient. The swimming strokes must 
be carried out with a characteristic frequency 
tJ^ to exploit rotational diffusion. Hence: 

6R D6R 
which gives, from Eq. (41): 



jjmigr 



DaSR 



Combining with Eq. (39), we obtain however: 

i?c D IR 

R ~ U'^^arR ^ z a' 



(42) 



Unless R ^ a and 5R/R is not too small (this 
would require a continuous swimmer undergoing 
large amplitude deformations) our device will 
have to swim a lot, before migration is able to 
overcome diffusion. 



9 Conclusion 

We have analyzed the behavior of a device that 

can swim by extracting energy from the gra- 
dients in an external shear flow. Adoption of a 
simple model, such as the; three-sphere swimmer 
of [2] , has allowed to identify optimal swimming 
strategies, both in infinite and wall bounded do- 
mains. 

In order to migrate across a shear fiow in an 
infinite domain, the microswimmer must main- 
tain on the average a configuration that is fore- 
aft asymmetric along the fiow. In order to mi- 
grate away from (towards) a wall perpendicu- 
lar to the direction of the shear gradient, the 
microswimmer must maintain on the average 
an elongated configuration along the stretching 
(contracting) direction of the external strain. 
Of these configurations, only the one with long 
axis along the stretching direction of the exter- 
nal strain, could be attained without an internal 
control system. In order for the extraction from 
the flow to take place, the microswimmer must 
maintain on the average by an elongated shape 
in a direction between — 7r/4 and ■n/A with re- 
spect to the flow. The energy extracted from 
the flow that is not dissipated by internal fric- 
tion, could be converted into swimming strokes, 
(and thus returned to the external flow), or in 
alternative, at least in principle, be stored in the 
swimmer, in the form of some potential energy. 

All the configuration described above could 
be obtained by a system of brakes controlling 
the stretching and contraction of the trimer arms 
in response to the external fiow. Inclusion of an 
elastic component in the dynamics may lead to 
higher swimming efficiency: we have shown that, 
in the case of an ideal trimer dynamics, propul- 
sion requires a deformation component for en- 
ergy extraction whose amplitude scales quadrat- 
ically in the amplitude of the swimming stroke; 
internal dissipation would cause this scaling to 
become linear. 

In the absence of an internal control sys- 
tem, a trimer, with dissipative springs between 
the beads, would be characterized, in a viscous 
shear flow, by the same orientation pattern as a 
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tank-treading vesicle [17] or a microcapsule [19]. 
The trimer would maintain, on the average, an 
elongated configuration, aligned with the flow 
in the case of infinite friction, and aligned with 
the stretching direction of the strain in the zero 
friction case. 

The natural scale for the migration velocity 
of a microswimmer in an external flow is the 

external velocity difference aR across its body, 
where a is the shear strength [see Eq. (1)] and 
R is the swimmer size. The swimming velocity 
jjmigr q£ ^j^g gwimmer in free space is going to be 
very small C/™97(ai?) = 0{a5R/R^), where a 
is the size of the moving parts (the beads) and 
5R is the amplitude of the swimming strokes. 
The correction by presence of a wall at distance 
L from the swimmer is going to be smaller by 
an additional factor {R/L)"^. 

Although small, the swimming velocity that 

can be achieved appear to dominate Brownian 
diffusion as long as the swimmer size remains 
above the micron threshold. Taking e.g. R = 
10/im. a = O.IR and SR = 0.3R, with a shear 
strength a = ls~^, we find a migration veloc- 
ity in free space C/^f^ ~ IQ-^cm/s [Eq. (14)] 
and a Brownian diffusivity D = KT/ {Girfxa) ~ 
10~^cm^/s, corresponding to migration domi- 
nating diffusion at distances above Rc ^ O.IR 
[Eq. (39)]. A swimmer that is ten times smaller, 
would be characterized in the same flow by {J™*^'" 
10~^cm/s for a diffusivity D ~ 10~*cm^/s and 
a crossover length Rc ~ lO'^R. Of course, the de- 
terministic theory breaks down in this regime, 
as Pe ~ 100 and the contribution to trimer de- 
formation from thermal noise becomes 0{1) [see 
Eqs. (34) and (38)]. In this regime, however, ro- 
tational diffusivity of the trimer could be ex- 
ploited to allow propulsion with use of just one 
degree of freedom (but with an internal motor: 
no more passive propulsion), which provides an- 
other example of swimming strategy to which 
the limitations of the scallop theorem do not 
apply. 

Coming back to the issue of swimming effi- 
ciency, as shown in [16], it appears that a swim- 
mer with a continuous body could achieve much 
better results than the ones described in the 
present paper, namely, 0{SR/R) rather than 
0{a6R/ R'^) efficiency. Notice that this is also 
better than the 0{{5R/ R)"^) result for a similar 
continuous swimmer in a quiescent fluid [29]. 
For SR/R ^ 1, the migration velocity JJ'^'^sr 
~ aR would have the necessary magnitude, to 



produce phenomena, analogous to the Fahraeus- 
Lindqvist effect in small blood vessels [38] . 

Appendix. The image field 

The image field must obey the equations of low 
Reynolds number hydrodynamics: 

VP = mV^U^ V • = 0, (Al) 

that is the Stokes equation plus incompressibil- 
ity, where P is the pressure. We can express U''^ 
in terms of scalar and vector potentials: 

= V# -I- V X A, (A2) 

where 

V^^ = and V • A = 0. (A3) 

The first of (A3) is a consequence of incompress- 
ibility; the second is a gauge condition. From 
here, the vorticity equation V x V^U'^ = 0, 
which descends from Eq. (Al), can be written 
in the form 

V^V^A = 0. (A4) 

Fourier transforming with respect to X2,z, the 
gauge condition can be written in the form 

A2k = -^A3k + f Alk, (A5) 

where the prime indicates derivative with re- 
spect to X2. The vorticity equation (A4), in- 
stead, becomes {d^^ — A:^)-^Ak = 0, whose gen- 
eral solution reads, imposing finiteness at Xi — )• 

— oo: 

Ak(Xi) = Ak (Xi - L) exp(A:(Xi - L)) 

akCxp(fc(Xi -i)). (A6) 

The second contribution to right hand side of 
Eq. (A6) is a pure gauge term that does not 
affect U^, and will be disregarded. The first of 
Eqn. (A3), instead, gives for <P: 

<Z'k(Xi) = <Pk exp(fc(Xi - L)). (A7) 

Using Eqs. (A2) and (A5), the expression for the 
velocity correction becomes, in terms of Fourier 
components: 

fC//k = C^3k + fe^lk + ^L, 

< i/4 = ifca^ik - ^3k + i^2'/>k, (A8) 
[UL = -It^sk + i^'i'k - ife^ik + iks^k, 
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and, imposing the boundary condition U^(L|0, Fi)tsined from Eq. (4): 
-Uk(L|0,Fi): 



Ik 



U2k = A3- ifca^k, 



^f/3k = f A3k-|fAik-iA;3^k, 



where Uk = Uk(i|0, Fi) and use has been made 
of Eqs. (A6-A7). Solution of this system gives: 



<?k = 
iak 
ilk 



_ i[-fc3fc2?73k+2ifc2fcr/ik + fc^!72k] 



2P" 

i(fc2C/3k-fc3C^2k) 

2k 



(A9) 



Substitution of Eqs. (A9) together with Eqs. 
(A6-A7) into Eq. (A8), gives, at Xi = 0: 

t^iyO|0,Fi) = -[(l + fcL)[/ik 
+ifc2-LC/2k + iksLUsi^] exp(-A;L), 

where Uk = Uk(-i^|0, Fi). In order to obtain 
AU™^^ , we thus have to antitransform 

5x,?7/k(0|X,Fi)|x=o = [(1 + kL)U[^ 
+ik2LU^k + ^k3LU3k]eM-kL), (AlO) 

where UJ^ = d\Jk/dL. Since the trimcr motion 
is confined in the X1X2 plane, we do not need to 
calculate 8x3 Uly^. We thus get, antitransforming 



Eq. (AlO) andax2C//k(0|X,Fi)|x=o 
(0|0,Fi): 

5xiC//(0|X,Fi)|x=o 

d^k 



(2^)2 y d^ilexp(-ik.Y^ 
X [(1 + kL)Ui + ik2LU^ + ik^LU'^], 



-ifc2C//k 

kL) 
(All) 



and 



a;C2C>/(0|X,Fi)|x=0 



= / ^ / d^Fi exp(-ik • Yx - kL) 

- kiLU2 - k3LU3]k2, (A12) 



(27r) 
x[i{l + kL)Ui 



where U = U(Y|0,Fi), Y = (L, Y^) and U' = 
dXj/dL. The stokeslet field at the wall U is ob- 



3a 



4al(L2 + y2)i/2 



U3 



_ 3a f 



Fv. 



+ 



4o-l(L2+x2)i/2 

3a Fa [^f 11+^2^12] 
4a (L2+y;2)3/2 



+ ^2^12] 

(L2+r2)3/2 

Y2[LFii + l2i^i2] 
(L2 + y2)3/2 



}• 



The integrals in Eqs. (A11-A12) arc carried out 
in polar coordinates with the help of MAPLE, 
and the result is Eq. (16). 
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